Endohedral terthiophene in zigzag carbon nanotubes: Density functional calculations 
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The inclusion and encapsulation of terthiophene (T3) molecules inside zigzag single- walled carbon 
nanotubes (CNTs) is addressed by density functional calculations. We consider the T3 molecule 
inside six semiconducting CNTs with diameters ranging approximately from 8 to 13 A. Our results 
show that the T3 inclusion process is exothermic for CNTs with diameters larger than 9 A. The 
highest energy gain is found to be of 2 eV, decreasing as the CNT diameter increases. This notable 
effect of stabilization is attributed to the positively-charged CNT inner space, as induced by its 
curvature, which is able to accommodate the neutral T3 molecule. The band structure of the 
T3@CNT system shows that T3 preserves its electronic identity inside the CNTs, superimposing 
their molecular orbitals onto the empty CNT band structure without hybridization. Our results 
predict that the electronic states added by the T3 molecules would give rise to optical effects and 
nonradiative relaxation from excited states. 



I. INTRODUCTION 

Filling carbon nanotubes (CNTs) with chosen 
molecules which in the gas phase exhibits particular elec- 
tronic properties are becoming increasingly important 
owing to the possibility they may modulate CNT proper- 
ties. Theoretical investigations have shown that semicon- 
ducting CNTs can be amphoterically doped in a control- 
lable way by encapsulating electrophilic and nucleophilic 
organometallic molecules^. Other calculations have sug- 
gested that confinement of molecules inside CNTs can 
induce structural phase transitions not seen in the bulk 
phase. For instance, encapsulated water in an 8 A- 
diameter CNT would exhibit a tight hydrogen-bonding 
network 2 as well as acetylene in a 7 A-diameter CNT 
would suffer polymerization. 3 In addition, recent exper- 
iments have demonstrated that template-grown large- 
diameter CNTs filled with fluorescent nanoparticles do 
emit fluorescent light, allowing the measurement of filling 
rates and transport of such nanoparticles inside CNTs. 

Endohedral functionalization with molecules needs a 
procedure to open the nanotubes at a first stage and 
then fill it with a particular molecule. Chemical methods 
are widely used, and the most simple one is using nitric 
acid which is very selective in attacking the carbon pen- 
tagon at the tube ends. 5 Recently, successful methods to 
fill opened multi-walled CNTs with fullcrcncs 6 and oth- 
ers organic molecules^ have allowed endohedral inclusion 
in a soft low-temperature way. Organic guest molecules 
like oligo and polythiophenes are an interesting family of 
functional electron systems which have been extensively 
studied recently due to their use in optical devices pSaLifii 
In particular, oligothiophenes (Tn) have spectral prop- 
erties dependent on the number of thiophene units (n). 
They show spectral changes and efficient energy trans- 
fer processes when included as guests in one-dimensional 
supramolecular architectures, as in perhydrotriphenylene 
crystals (PHTP)»±i This characteristic can be exploited 
for energy conversion or in detectors of short-wavelength 
radiation. Recently, we have reported conformational de- 



tails of the PHTP:T3 inclusion compound^ 

In this work, we study the energetic and electronic 
properties of single-walled CNTs filled with terthiophene 
(T3) molecules as well as the T3 incorporation mecha- 
nisms by a open-ended CNT. Single-walled CNTs, which 
are typically 10 A in diameter, can be metallic or semi- 
conducting depending on the way that a graphene sheet 
is rolled up, which is characterized by the chiral indices 
(n, m). Metallic tubes occur if n — m is divisible by 3; oth- 
erwise, the tubes are semiconducting M> Non-chiral nan- 
otubes with indices (n, 0) and (n, n) are termed zigzag 
and armchair, respectively, which is related to the ar- 
rangement of the C atoms around the tube. We study 
the T3 molecule inside six zigzag CNTs, with diameters 
larger than 8 A. Our results show that the inclusion of 
the T3 molecule is an exothermic process for all CNTs 
under consideration and that the energies involved de- 
pend on the diameter of the nanotube. We also found 
that the highest occupied and lowest unoccupied molec- 
ular orbitals (HOMO and LUMO) of the free T3 molecule 
appear as sharp subbands when these molecules fill the 
tube, showing additional singularities in the density of 
states with similar energy difference than that of the free 
molecule. This result suggests that the T3 optical prop- 
erties could be preserved when this molecule fills a semi- 
conducting CNT. 



II. THEORETICAL METHOD 

The calculations were carried out in the framework 
of density functional theory (DFT)^4 with the gener- 
alized gradient approximation (GGA) to the exchange- 
correlation functional^ We use a basis set consisting of 
strictly localized numerical pseudoatomic orbitals, as im- 
plemented in the SIESTA code^ namely double-^ singly 
polarized basis set. Norm-conserving pseudopotentials^i 
in their separable form 18 are used to describe the 
electron-ion interaction. We consider the zigzag (n,Q) 
CNTs with n=10, 12, 13, 14, 15, and 16, which have 
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FIG. 1: (Color on line) Atomic geometry of the T3 molecule 
in the gas phase. The three thiophene units adopt a nonpla- 
nar configuration (lower figure) characterized by the dihedral 
angle S-C-C-S (0 D ). 



diameters of 8.0, 9.6, 10.2, 11.1, 11.4, and 12.7 A, re- 
spectively. Infinite CNTs are described within the su- 
percell approach with periodic boundary conditions ap- 
plied only along the nanotube axis. As the T3 molecule 
has a length of about 12.4 A, we choose a periodic CNT 
unit of 17.3 A in length. Thus, the distance between T3 
molecules inside a nanotube (T3@CNT) is about 5 A, 
ensuring a negligible interaction between them. We con- 
sider a lateral distance between tubes of about 8 A. We 
also study the T3 incorporation through an open-ended 
CNT. This was done by taking a finite CNT segment of 
17.3 A plus the T3 molecule which is located initially 
2 A apart from the tube end, along its axis. Then, the 
molecule is gradually moved into the tube. We consider 
two possible scenarios for the CNT entrance, with bare C 
atoms and those saturated by H atoms. The other CNT 
end is always saturated by hydrogen. The positions of 
all atoms in the nanotubes and the molecule were re- 
laxed using the conjugated gradient algorithm until the 
force components become smaller than 0.05 eV/A. Due 
to the large size of the supercells, we use the T point for 
the Brillouin zone (BZ) sampling. We check the band- 
structure calculations and the equilibrium geometries for 
the T3@(12,0) and T3@(14,0) systems considering four 
k points, according to the Monkhorst-Pack meshiiS Our 
results do not show any significant variation with respect 
to the T-point calculations. We also check the difference 
in energy between T3@(12,0) and T3@(14,0), consider- 
ing one and four points. We find a very small energy 
difference of 3 meV/atom, showing that the T point is 
sufficient to ensure the convergence of BZ sampling. 



III. RESULTS AND DISCUSSION 

We first discuss the structural and electronic prop- 
erties of the isolated molecule and the semiconducting 
CNTs. Gas phase T3 (Ci2H 6 S3) has been theoretically 
characterized by a S-C-C-S dihedral angle of 147.2°^ 
and experimentally by the Sq — > Si singlet transition 
of 3.05 eVi^i Figure n shows our results for the equilib- 
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FIG. 2: (Color on line) Energy profile during the axial inser- 
tion of the T3 molecule into the (13,0) CNT by both clean 
and hydrogenate entrances. The inset shows the T3 molecule 
entering inside the (13,0) CNT by the clean entrance. 



rium geometry; we have obtained a more planar struc- 
ture with a dihedral angle of 157.8°, whereas the T3 
HOMO-LUMO energy is found to be 2.33 eV. This dis- 
crepancy is a consequence of the DFT-GGA methodol- 
ogy which underestimates the gap energy. Concerning 
the carbon nanotubes, we find that (12,0) and (15,0) 
CNTs are small-gap semiconductors with a gap energy 
of 0.04 and 0.05 eV, respectively, in close agreement with 
experiments^ whereas (10,0), (13,0), (14,0) and (16,0) 
CNTs have a larger band gap of 0.85, 0.70, 0.60, and 
0.60 eV, respectively. 

Figure [5] shows the variation in the total energy when 
T3 enters into a (13,0) CNT by a clean and hydrogenated 
entrance. We find similar energy profiles in both cases in- 
dicating that electrostatic forces induced by the dangling 
bonds at the CNT termination are not determinant for 
the insertion process. We observe that the incorporation 
reaction is exothermic by 1.2 eV. Thus, it is energetically 
favorable that the molecule enters into an uncapped nan- 
otube than remains outside, exhibiting a barrierless in- 
sertion process. From the slope of the curves of Fig. [21 we 
can estimate the capillary forces involved in the molecule 
insertion, considering that F = AE/Ad. We find that 
the force which pushes T3 into the tube is about 0.15 
nN. In a recent work, Kim et al. studied the insertion of 
a C2H2 molecule into the (5,5) armchair CNT which has 
a diameter of about 7 A. They found a force of 0.22 nN, 
very close to our results even though the C2H2 molecule 
is smaller and has a linear geometry^ This suggests a 
general trend for the molecular insertion process in car- 
bon nanotubes. 

We investigate the energetic of zigzag CNTs of dif- 
ferent diameters filled with T3 molecules. In our su- 
percell approach, the molecules are aligned coaxially in- 
side the tube with distances separating them of about 



3 



3.0 rn — | — i — i — i — i — | — i — i — i — i — | — r 



2.5 



> 
u 2.0 h 



ta 1.5 - 

so 

g 1.0 



0.5 h 
0.0 



~1 — i — i — i — i — I — i — i — i — r 




(10,0) 



(15,0) 



(16,0) 



J I I I I I I I I I I I I I I I I I I I I I I I I l_ 



9 10 11 12 

Tube diameter (A) 



13 



FIG. 3: (Color on line) Energy gain as a function of the CNT 
diameters for the T3 molecule inside the (n,0) CNTs, with 
n = 10 — 16. The zero energy corresponds to T3 in the gas 
phase plus the empty CNT. The inset shows the unit cell of 
T3@(13,0). 



5 A. We find that the energy of the T3@CNT system 
is lower than the isolated constituents. Figure [31 shows 
the gain in energy of the T3@CNT systems with respect 
to the empty CNT plus the isolated T3 molecule. We 
find that T3@(12,0) has the highest energy gain of about 
2 eV, whereas increasing the CNT diameters, this en- 
ergy tends to decrease. This notable effect of stabiliza- 
tion is attributed to the difference in the relative charge 
between the inner and outer surfaces of the CNT. For 
small CNTs, the higher curvature produces a polariza- 
tion normal to the tube surface which tend to dislocate 
the electronic charge toward outside, as shown by Du- 
mitrica, et al., based on a rehybridization model and DFT 
calculations. a» Thus, a neutral molecule with a negative 
electronic cloud surrounding the positive nucleus would 
prefer to accommodate in the slightly positive CNT in- 
ner region than remain outside. Therefore, nanocapil- 
larity forces acting on the T3 molecules would originate 
in the normal polarization between the inner and outer 
CNT surfaces. For larger-diameter CNTs the polariza- 
tion tends to disappear because the inner and outer cur- 
vatures tend to equalize. It is worthy to note that the 
gain in energy for the neutral T3 molecule inside CNTs 
can be as high as 1.5 eV for a small diameter variation of 
about 3 A, which is the difference in diameters between 
(12,0) and (16,0) CNTs. 

Our results for the T3 molecule inside the small (10,0) 
CNT of 8 A in diameter show a gain in energy of 
1.1 eV with respect to the isolated components. How- 
ever, the CNT shows a strong oval deformation whereas 
the molecule becomes completely flat, resulting in a tight 
encapsulation, as shown in Fig. Ufa). Although the (10,0) 
CNT has a stronger normal polarization owing to its en- 
hancing curvature, the T3 inclusion process appears to 




FIG. 4: (Color on line) Atomic geometry of the T3 
molecule inside carbon nanotubes of different diameters, (a) 
T3@(10,0), (b) T3@(12,0), and (c) T3@(13,0). The black dot 
indicates the nanotube centers. 



be unrealistic due to steric effects. The above result sug- 
gests that benzenelike molecules would not enter into un- 
capped zigzag CNTs with diameters closer to 8 A. 

Table [I] compares the geometry and electronic proper- 
ties of T3 in the gas phase and inside the CNTs. We 
find a larger molecular deformation in T3@(12,0), reduc- 
ing the diahedral angle in 14° with respect to the free 
molecules. This tendency reverses when the CNT diam- 
eter increases. Figures Etb) and^c) show an axial view 
of the T3@(12,0) and T3@(13,0) systems. We observe 
that the T3 molecule inside the (12,0) CNT is shifted 
about 0.5 A toward the wall from the CNT axis with the 
S atoms facing the wall, whereas inside the (13,0) CNT 
the molecule remains at the CNT axis, increasing the 
dihedral angle. This result somewhat confirms the posi- 
tive character of the tube inner wall which in the (12,0) 
CNT induces a sufficiently high normal polarization that 
moves the molecule toward the wall, decreasing the dihe- 
dral angle. Table U also shows that both the CNT band 
gap and the T3 HOMO-LUMO energy remain essentially 
unchanged when the molecule is inside the CNTs. 

Figure shows the density of states (DOS) of the 
empty CNTs (13,0), (14,0), and (16,0) (upper panels), 
which are compared with the DOS of the same CNTs 
filled with T3 molecules (lower panels). The DOS of 
filled CNTs clearly shows the appearance of both HOMO 
and LUMO states of T3 at essentially the same energy, 
about Ep — 0.44 eV and Ep + 1.86 eV, respectively. 
The difference in energy between these states is found 



TABLE I: Structural and electronic parameters for the T3 
molecule and the T3QCNT system. 9o is the dihedral angle, 
and Ehl is the HOMO-LUMO energy of the molecule. Eg is 
the T3@CNT gap energy. 



System 


6> D (deg) 


£hl (eV) 


E G (eV) 


T3 


157.8 


2.33 




T3@(12,0) 


143.8 


2.30 


0.05 


T3@(13,0) 


157.3 


2.29 


0.70 


T3@(14,0) 


155.0 


2.34 


0.60 


T3@(15,0) 


167.8 


2.34 


0.05 


T3@(16,0) 


161.8 


2.27 


0.60 
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FIG. 5: Density of states (DOS) of empty CNTs (13,0), (14,0), 
and (16,0) (upper panels) and DOS of the same CNTs filled 
with T3 molecules (lower panel). The arrows indicate the 
HOMO-LUMO energy width of T3. 



to be of 2.30 eV, very close to that calculated for the 
free molecule, 2.33 eV. Other molecular states can be 
observed between E F - 2.0 eV and E F - 1.5 eV. This 
suggests a minor influence of the inner CNT electrostatic 
potential in the T3 molecular states. Therefore, we may 
speculate that the optical properties of the T3 molecule 
would be maintained inside the tube. Nevertheless, the 
addition of molecular states to the band structure of the 
CNT itself produces a system with more complex optical 
characteristics, where we may predict fluorescence from 
the T3 molecule as well as fluorescence due to transi- 
tions between the T3 HOMO band to the sharpest Van 
Hove singularity in the conduction band. In addition, 



highly probable energy transfer processes between olig- 
othiophcnc molecules and the CNT may occur. In this 
case, a plausible mechanism would allow absorption in 
the near UV (Sq — » Si of T3), with energy transfer 
from T3 (acting as donor) to the CNT (acting as accep- 
tor). For instance, down-conversion processes between 
Van Hove singularities would give fluorescent emission 
from up to the En IR transition, in the case of the 
T3@(13,0) system. It is important to remember that the 
calculated energy transitions were obtained within the 
DFT-GGA approach which underestimates their values. 
Thus, the optical properties extracted from the above 
results must be taken as a qualitative picture. 



IV. SUMMARY AND CONCLUSIONS 

A notable effect of stabilization is produced by the en- 
dohedral inclusion of the oligomer through open-ended 
zigzag single-walled CNTs, which is the order of magni- 
tude of a chemical bond. This effect is observed with 
similar magnitude in both bare and hydrogenated CNT 
entrances. Earlier studies of several neutral fullerenes 
predict that the electrostatic potential inside the cage 
has a positive charge character and consequently would 
be suitable to accommodate anionic or neutral species^ 
More recent calculations confirms this finding in carbon 
nanotubes^ showing that the curvature of the CNT 
walls produces a charge displacement toward the out- 
side, leaving a positive character inside. Thus, the T3 
inclusion or encapsulation in a CNT is driven in a sig- 
nificant extent by nanocapillarity forces induced by the 
interaction between the slightly positive CNT cavity and 
the electronic cloud of the molecule, giving an exother- 
mic process. Thus, the drop in the energy gain for the 
T3 inclusion in larger-diameter CNTs can be explained 
in terms of their decreasing curvature. 

The band structures of T3@CNTs show that the T3 
molecule preserves its electronic properties, superimpos- 
ing the HOMO-LUMO states onto the original CNT band 
structure. Semiconducting single-walled CNTs are opti- 
cally active systems that are fluorescent when encased in 
cylindrical micelles with a striking influence on the diam- 
eter and chiral angle^*2^ showing strong absorption and 
fluorescence due to their sharp Van Hove singularities. 
Our results predict that the electronic states added by 
the T3 molecules to the CNT band structure would give 
rise to optical effects in the radiative relaxation from the 
excited states. In addition, energy-transfer processes be- 
tween T3 and the CNT suggest possible applications in 
optical nanodevices, for instance, in UV down-converters. 
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